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ABSTRACT 

Although abundance gradients in the Milky Way Galaxy certainly exist, de- 
tails remain uncertain, particularly in the inner Galaxy, where stars and H II 
regions in the Galactic plane are obscured optically. In this paper we revisit two 
previously studied, inner Galaxy H II regions: G333.6— 0.2 and W43. We ob- 
served three new positions in G333.6— 0.2 with the Kuiper Airborne Observatory 
and reobserved the central position with the Infrared Space Observatory’s Long 
Wavelength Spectrometer in far-infrared lines of S ++ , N ++ , N + , and 0 ++ . We 
also added the N + lines at 122 and 205 /rm to the suite of lines measured in W43 
by Simpson et al. (1995). The measured electron densities range from ~ 40 to 
over 4000 cm -3 in a single H II region, indicating that abundance analyses must 
consider density variations, since the critical densities of the observed lines range 
from 40 to 9000 cm -3 . We propose a method to handle density variations and 
make new estimates of the S/H and N/H abundance ratios. We find that our sul- 
fur abundance estimates for G333.6— 0.2 and W43 agree with the S/H abundance 
ratios expected for the gradient previously reported by Simpson et al., with the 
S/H values revised to be smaller owing to changes in collisional excitation cross 
sections. The estimated N/H, S/H, and N/S ratios are the most reliable because 
of their small corrections for unseen ionization states (< 10%). The estimated 
N/S ratios for the two sources are smaller than what would be calculated from the 
N/H and S/H ratios in our previous paper. If all low excitation H II regions had 
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similar changes to their N/S ratios as a result of adding measurements of N + to 
previous measurements of N ++ , there would be no or only a very small gradient 
in N/S. This is interesting because nitrogen is considered to be a secondary ele- 
ment and sulfur is a primary element in galactic chemical evolution calculations. 

We compute models of the two H II regions to estimate corrections for the other 
unseen ionization states. We find, with large uncertainties, that oxygen does not 
have a high abundance, with the result that the N/O ratio is as high ( 0.35) 

as previously reported. The reasons for the uncertainty in the ionization correc- 
tions for oxygen are both the non-uniqueness of the H II region models and the 
sensitivity of these models to different input atomic data and stellar atmosphere 
models. We discuss these predictions and conclude that only a few of the latest 
models adequately reproduce H II region observations, including the well-known, 
relatively-large observed Ne ++ /0 ++ ratios in low- and moderate-excitation H II 
regions. 

Subject headings: Galaxy: abundances — H II regions — ISM: individual(G333.6— 0.2, 
W43) 


1. Introduction 

The presence of radial abundance gradients in the plane of the Milky Way is now an 
established fact and is seen in both stars and gaseous nebulae (see, e.g., Henry & Worthey 
1999; Rolleston et al. 2000; and references in both papers). The cause is generally thought 
to be due to the formation history and subsequent evolution of the Galaxy; thus the observed 
gradients are a major tool for understanding this history (see, e.g., Hou, Prantzos, & Boissier 
2000; Chiappini, Matteucci, & Romano 2001; Chiappini, Romano, & Matteucci 2003, and 
references therein). The basic idea is that the inner Galaxy formed before the outer Galaxy 
and the higher molecular gas density in the inner Galaxy produces a higher star formation 
rate. The result is a greater return to the interstellar medium (ISM) in the inner Galaxy 
of both “primary” alpha elements from massive star supernovae and “secondary” elements 
like nitrogen. Secondary nitrogen is produced by CNO burning of already existing carbon 
and oxygen in intermediate-mass stars and subsequently returned to the ISM through mass 
loss. However, both the wide range of uncertain input parameters to chemical evolution 
models (Pagel 2001) and the uncertain details of the abundance variation of each element, 
primary and secondary, contribute to our less-than-complete understanding of the formation 
and evolution of the Milky Way. In this paper we try to improve our knowledge of the 
abundances in inner Galaxy H II regions and also test the methods that have been used to 
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date to determine abundance gradients from low-excitation, optically obscured H II regions. 

H II regions in the inner Galaxy (galactocentric radius Rg < 6 kpc, where Rq = 8 kpc) 
are not accessible to optical observers because of interstellar extinction. Consequently, the 
abundances must be determined from far-infrared (FIR) lines, such as those of N + at 122 
and 205 pm, N ++ at 57 //m, 0 ++ at 52 and 88 pm, Ne + at 12.8 pm, Ne ++ at 15.5 and 36 
pm, S ++ at 18.7 and 33.5 pm, and S 3+ at 10.5 pm. Large surveys include those of Simpson 
et al. (1995a, hereafter SCREH) and Afflerbach et al. (1997), who observed the FIR lines 
(wavelengths > 17 pm) with the Kuiper Airborne Observatory (KAO) in 18 inner Galaxy 
H II regions, and Martin-Hernandez et al. (2002a), who observed 13 H II regions with the 
Infrared Space Observatory (ISO). The KAO observations included FIR measurements of 
the lines with wavelength > 18 pm, whereas the ISO observations included the mid-infrared 
lines, measured with ISO’s Short Wavelength Spectrometer, SWS, but not the [N II] 205 pm 
line. 


A major result of these studies is that inner Galaxy H II regions have low excitation 
compared to typical H II regions at larger Rg-, this is apparent in both heavy element ionic 
ratios such as 0 ++ /S' t " f (SCREH) and Ne ++ /Ne + (Simpson & Rubin 1990; Giveon et al. 
2002), and also in He + /H + measured from radio recombination lines (Churchwell et al. 1978; 
Thum, Mezger, & Pankonin 1980). It is still controversial whether the observed increase in 
excitation with increasing Rg is due entirely to heavy element opacity effects in the H II 
regions and stellar atmospheres, or whether there is an additional gradient in the maximum 
stellar effective temperature, T e ff, of the exciting stars as well (e.g., Giveon et al. 2002; 
Martin-Hernandez et al. 2002b; Smith, Norris, & Crowther 2002; Morisset et al. 2003). 

Two of the abundance gradients that are most needed for galactic chemical evolution 
studies are the O/H and N/H ratios, representing primary and secondary elements, respec- 
tively. In optically obscured H II regions the only ionization state of oxygen available for 
study is 0 ++ , which means that corrections are required for the unseen 0 + ions. The sit- 
uation is better for nitrogen, as there are [N II] lines as well as [N III] lines in the FIR 
wavelength range. Unfortunately, the [N II] 122 pm line is difficult to measure for airborne 
observations because of atmospheric absorption, and the [N II] 205 pm line is difficult because 
of poor detector sensitivity. In addition, the [N II] lines are difficult to interpret because 
they are collisionally de-excited at much lower densities than the lines from the other ions. 
Consequently, the papers to date have used only the [N III] and [O III] lines to estimate N/O 
and the secondary/primary abundance ratio in optically obscured H II regions (SCREH; 
Afflerbach et al. 1997; Martin-Hernandez et al. 2002a). This is particularly troubling since 
a majority of the nitrogen and an even larger fraction of the oxygen are singly ionized in the 
high metallicity H II regions of the inner Galaxy (note that in the outer Galaxy, the lower 
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metallicity H II regions can be much more highly excited and thus the N ++ /0 ++ ratio is 
often a good approximation to N/O, Rudolph et al. 1997). On the other hand, a majority 
of the primary element sulfur is doubly ionized in H II regions. Here we devise a method of 
analyzing the density-sensitive [N II] lines to more accurately estimate nitrogen abundances 
with respect to sulfur without the need for large corrections for unseen ionization states. 
Using the observed N ++ /N + ratio and calculated H II region models, we then estimate the 
ionization correction factors (ICFs) for oxygen. 

Two of the most luminous H II regions in the Galaxy are G333.6-0.2 at Rq ~ 5.5 kpc 
(Wilson et al. 1970) or 5.6 kpc (Caswell & Haynes 1987) and W43 (G30.8— 0.0) at Rq ~ 4.2 
kpc (Reifenstein et al. 1970) for Rq = 8.0 kpc. The radio fluxes from these objects require 
ionizing stellar luminosities of 10°° and 2.9 x 10 50 photons s -1 , respectively (see Rubin 
1968b). Both H II regions are located in the Scutum-Crux spiral arm (see Vallee 2002) in 
the Milky Way molecular ring. Both are excited by clusters of massive stars, observed at 
near-infrared (NIR) wavelengths. The cluster exciting G333.6— 0.2 contains stars that are 
very young, embedded in dust (Blum et al. 2002), in keeping with the youthful appearance 
of this extremely dense, compact H II region. On the other hand, the cluster exciting W43 
must be older and more evolved because it contains a Wolf-Rayet star in addition to O giants 
or supergiants (Blum, Damineli, & Conti 1999; Cotera & Simpson 1997), although there are 
still numerous stars forming in the associated molecular clouds (Motte, Schilke, & Lis 2003). 
This is also consistent with the more mature appearance of W43, where the ionized gas is 
well separated from the exciting star cluster. 

In this study, we report measurements of the [N II] lines in the three positions observed in 
W43 by SCREH and numerous lines in three additional positions in G333.6— 0.2 north of the 
center position observed by Colgan et al. (1993) with the KAO; we also report measurements 
of the central position in G333.6— 0.2 with ISO’s Long Wavelength Spectrometer (LWS). In 
Section 2 we describe the observations, in Section 3 we discuss analysis techniques that are 
used to determine abundances in both low- and high-density H II regions, in Section 4 we 
model the two H II regions to determine the ICFs needed for those elements without measured 
abundances of the singly ionized atoms, and in Section 5 we discuss models in general and 
give our conclusions. By measuring both N + and N ++ , we hope to avoid the need for large 
(and uncertain) ICFs for nitrogen and thus improve the reliability of the secondary/primary 
abundance ratio measurement in the Milky Way molecular-ring H II regions. 
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2. Observations 

Both G333.6— 0.2 and W43 are extended compared to our 45-60" KAO beam and 80" 
LWS beam, which means that integrating over a constant density model to calculate ICFs is 
not adequate. Instead, we measured the fluxes at a number of positions at varying distances 
from the exciting stars to compare with the radial distributions of line flux computed from 
models to determine the source excitation and abundances. 


2.1. KAO Observations 

The observations were made with the Cryogenic Grating Spectrometer (CGS, Erickson 
et al. 1985) on NASA’s 0.91 m-KAO telescope. Data acquisition and reduction techniques 
are described by Colgan et al. (1993) and SCREH. For the present observations, Saturn was 
used for calibrating the absolute fluxes and determining the relative detector responses. The 
observed line and continuum fluxes for all positions are given in Table 1. 

G333.6— 0.2 was measured at its center and in three offset positions. The data from 
the center position, C (RA 16 h 22 m 9!3, Dec — 50°6'5", both J2000), were taken with a 45" 
beam and were previously reported by Colgan et al. (1993). The three North positions (Nl, 
N2, and N3 at 1', 2', and 3' north of G333.6— 0.2C) were measured on 1993 May 11, 13, 
and 20, flying from Christchurch, New Zealand. For these positions, a 60" beam was used 
for all lines to provide good ratios to the [N II] 205 pm line, at whose wavelength the KAO 
has a diffraction size of 1.22X/D = 57". The chopper angle was approximately east-west 
(position angle ~ 90°) and the throw was ~ 8.5'. Because the [N II] 122 pm line is on the 
edge of a strong telluric H 2 0 line, the observations of positions Nl, N2, and the [N II] 205 
pm line in Position N3 were taken at or near 43,000 feet altitude, instead of the more usual 
39,000-41,000 feet, where the other lines in Position N3 were measured. Corrections for 
telluric H 2 0 absorption were made for both source and calibrator — we estimate the zenith 
H 2 0 overburden to range from 4.4 to 5.2 pm at 43,000 feet on May 11 and 13, and 3.9 pm 
at 41,000 feet on May 20. 

The three positions of W43 that had previously been studied by SCREH (W43N at 
RA 18 h 47 m 36?8, Dec -1°55'29"; W43C at RA 18 h 47 m 36 s 0, Dec -1°56'39"; W43S at RA 
18 h 47 m 38®5, Dec — 1°57'34", all J2000) in the [N III], [O III], and [S III] lines were measured 
in the [N II] lines on 1992 July 1, flying from NASA Ames. The 122 pm line was measured 
in both 45" and 60" beams and the 205 pm line was measured in a 60" beam. The chop- 
per angle was approximately east-west, the throw was ~ 6', and the zenith telluric H 2 0 
overburden was ~ 5 pm at 43,000 feet and ~ 2 pm at 45,000 feet. We estimate that there 
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are extra measurement uncertainties of order 20% in addition to the statistical errors given 
in Table 1 for the W43 [N II] 122 p m lines. These could result from dumpiness, pointing 
differences, or the effects of poor correction for deep atmospheric absorption. The analyses 
of the various measurements produce different abundance measures, as will be tabulated in 
the next sections. 


2.2. ISO Observations 

G333.6— 0.2 was measured with the ISO Long Wavelength Spectrometer in LWS01 mode 
on 1997 March 20 (TDT No. 49001002, PI R. Rubin), exposure time 1054 s. The line 
strengths for the [O III] 52 and 88 /. im and [N III] 57 gm lines are given in Table 1 (the [O I] 
63 gra line is also strong, but the [N II] 122 pm line was not detected owing to ISO's low 
spectral resolution in LWS01 mode with the LW detectors). 


2.3. Extinction Corrections 

The measured line fluxes were corrected for extinction using the infrared extinction 
function of Li & Draine (2001). This function has a very similar wavelength dependence 
to the extinction function used by SCREH when normalized to Li & Draine’s extinction 
maximum at ~ 9.5 pm. Thus the FIR optical depths are characterized as being proportional 
to the maximum optical depth in the 9.5 to 10 pm region, t 9 . 7 . Based on new information in 
the literature, we have updated these values from SCREH, with W43 having more extinction 
and G333.6— 0.2 having less. 

For three stars in the central cluster of W43, Blum et al. (1999) estimated that Ax — 
3.55, 3.52, and 3.63. Using the extinction law of Li & Draine (2001), we find that Ax = 3.55 
corresponds to r 9 . 7 = 2.74. Thus we increase the values of r 9 . 7 used by SCREH by 37% to 
r 9 7 = 2.74, 2.74, and 3.63 for W43N, W43C, and W43S, respectively. 

For G333.6— 0.2, Fujiyoshi et al. (2001) inferred an average value of r 9 . 7 = 1.5 from their 
mid-IR spectropolarization measurements. Consequently, we use r 9 . 7 = 1.5 for the center 
position instead of SCREH’s r 9 . 7 = 2.0. (This value corresponds to T Hq = 14 mag; a lower 
extinction yet for the center of G333.6— 0.2 might be predicted, since the nebula was detected 
in Ha by Churms et al. 1974.) We use r 9 . 7 = 1.5 for the North positions as well, since there 
are no significant variations in the integrated H I optical depths measured by Forster et al. 
(1987) at the C, Nl, and N2 positions. 

Following SCREH, an extinction uncertainty estimate of 20% (la) is propagated through 
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all subsequent data analysis errors. 


3. Analytical Approach 

Because recombination, bremsstrahlung, and collisional excitation rates are all propor- 
tional to the electron density squared, N%, it has long been the practice to treat all abundance 
analyses in H II regions and planetary nebulae (PNe) as though the surface brightness is, in 
fact, proportional to the emission measure, EM = f Nidi , where l is the line-of-sight path 
through the nebula. This is valid so long as N e is substantially less than the critical density, 
Ncrit, defined as the density at which the upper level of a transition is depopulated equally by 
radiative transitions and by collisions (e.g., Osterbrock 1989; Rubin 1989). Corrections for 
collisional de-excitation can be made if N e is known and approximately constant. However, 
constant densities are unlikely to occur in real nebulae. The effect of density variations is 
that the kernels in the integral over the nebular volume have different weighting functions 
for ions with different amounts of collisional de-excitation. Thus the flux ratios of such ions 
depend on these variable weighting functions (see, e.g., Rubin 1989). Here we consider the 
calculation of ratios of column densities of ions whose emission is proportional to N e , as well 
as to Nl . The former is appropriate when N e Ncrit- These two methods give identical 
results for ion ratios when the density is constant over the whole observed volume. We will 
demonstrate the effects of density variations via a simple example and use the results to 
estimate how to handle the transition case where N e ~ Ncrit- 


3.1. Low Density H II Regions 

For a given ion, a “low density” H II region is one where N e <C N^u and a “high 
density” H II region is one where N e Ncriu thus any particular H II region can have both 
low density and high density volumes, depending on the value of N crit for the transition 
being studied. Rubin (1989) tabulated values of N crit for a large number of optical and IR 
transitions and Martfn-Hernandez et al. (2002a). tabulated values for many IR transitions; 
these values, of course, need to be updated as the atomic data are revised. The values of 
Ncrit have a strong dependence on wavelength (owing to the dependence on wavelength of the 
transition probabilities), ranging from ~ 40 cm -3 for the [N II] 205 pm line to > 10 8 cm -3 
for some ultraviolet transitions. Because of their high Ncrit compared to typical densities 
observed in H II regions and PNe, most optical and near-infrared (NIR) forbidden line and 
recombination line transitions are analyzed in the low density limit. 
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In the optically thin case (no dependence on line profile function), the line emissivity, 
j, used in the radiative transfer equation can be written as 


j 


hv ^ N upper 

~a at 

47 r Ns 


(1) 


where A is the Einstein transition probability, Ni is the density of atoms in the particular ion- 
ization state, N upper /Ni is the fraction of ions in the upper level of the transition, and hv is the 
energy of the transition. Then the observed line flux (ignoring extinction and self-absorption) 
is F = f jdldQ where dldQ is the volume element of the H II region in the telescope beam, 
Q. We define the normalized volume emissivity e = 4nj/NiN e (Rubin et al. 1994; SCREH) 
and the average ^-weighted ionic density as < N{ >= f NiN e dldQ/ f N e dldVt. For the low 
density regime where N e N cr n, e has very little dependence on N e , and the line surface 
brightnesses (or fluxes) are proportional to f N^dl (or f N^dldQ). 


Thus for this case we can write 


F = 


/ 


—NiN e dldQ 

47T 


e 

4w 


< Ni > 


N e dldQ 


( 2 ) 


and 


Fl _ ^1 < > 

F2 62 < Ni 2 > ’ 


( 3 ) 


so the ionic abundance ratio is 

^ ^ __ -^1 / e i 

< N i2 > ~ F 2 /e 2 - 

As N e approaches N cr n, collisional de-excitation of the line’s upper level becomes more and 
more important, and the ratio of 2 lines from the same ionic species Ni with different N cr n 
can be used to estimate the electron density (e.g., Rubin 1989; Fig. la). 


3.2. High Density H II Regions 


At high densities, the fractional population of the upper level approaches that of local 
thermodynamic equilibrium (LTE) and is a function of electron temperature, T e , only. In 
this regime, the line surface brightness is proportional to the column density, CD = f Nidi , 
the number of ions in the line of sight: 


F = 


w 

4tt 


NidldQ rs 



Nidi, 


( 5 ) 


where the quantity N e e is the emissivity per ion and has very little density dependence so 
long as there are no low density regions along the line of sight (no large density variations). 
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Then for identical beam areas f2, the flux ratio is 

F\ (1 %e l )CD, 

F 2 (N k c 2 )CD 2 ’ ' 

and the column density ratio is 

£A = F,/W, fl ) 

CD 2 F 2 /(N„e 2 y 1 ' 

In fact, for abundance ratio calculations, one really wants to measure the actual number of 
ions in the line of sight, and not the number of ions weighted by some other, strongly varying 
parameter such as the electron density (as in equation (4)). Thus it is desirable to estimate 
column densities for gaseous nebulae whenever feasible. 


3.3. Density Variations: An Illustrative Model 

The presence of density variations has a number of effects on the abundances that one 
infers from FIR line measurements. We demonstrate these effects by the use of simple slab 
models, with an exponentially decreasing density N e = Noe~ x from a location in the interior 
of the slab x 0 to an outer location x mai . This density function was chosen because it gives 
similar weighting to both high and low density regimes (other types of models, e.g., Rubin 
1989, could be used to illustrate the same effects). For simplicity, the example assumes no 
ionization stratification, iV* = N e , and T e = 7000 K; thus the only effects on the calculated 
volume emissivities are those of the variation in density. For all models, a very large value 
of Zmax was chosen such that N e = 0.01 cm -3 . Each model had a different value of xo and 
A r 0 , ranging from N 0 = 1 cm -3 to N 0 = 10 6 cm -3 . The intensities emitted in each line were 
calculated from 

fXmax € ( r \ 

/= / -p- NiNedx . (8) 

Figure la show's how the electron density, N est , can be estimated using calculated volume 
emissivity ratios from the observed [S III] 19/33 /zm, [O III] 52/88 /zm, and [N II] 122/205 
/zm line ratios, wdth the [N II] ratio being sensitive to low' densities and the [S III] ratio 
sensitive only to much higher densities. Values of N est were computed from the ratios of the 
line intensities produced by the models; ratios of these estimated densities to the maximum 
density of each model, 7V 0 , are plotted in Figure lb. This figure show's that in the presence 
of density variations, even w'hen the ions coincide exactly in space, -one does not infer the 
same density from the three ratios: N est ( N II) < N est ( O III) < N est ( S III). From this, one 
must conclude that the volumes containing the observed ions are not co-extensive (as is the 
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case for, e.g., ionization stratification) if the estimated densities are in the range of density 
sensitivity and do not follow these inequalities (subject to the observational and atomic data 
uncertainties). 


To understand the effects of density variations on estimated abundances, we calculate 
the estimated numbers of ions in the volume from the model line intensities and N es t and 
compare these estimates with the true numbers of ions found by integrating N e or TV 2 over 
the model line of sight. These comparisons we call “abundance measures”. Numerically, 
we define the abundance measures, AM, as either the calculated line intensity / divided by 
the normalized volume emissivity e for low density analysis or / divided by the emissivity 
per ion (N est e) for high density analysis. These quantities are normalized by dividing by 
the appropriate integrals of density squared or density, respectively. That is, the abundance 
measure for analysis that assumes that the line intensities are proportional to the emission 

measure is 

I 


AM EM — — 
e 


NiNpdx 


r Xq 


-1 


( 9 ) 


and the abundance measure for analysis that assumes that the line intensities are propor- 
tional to the column density is 


AMcd — 


I 

( N est e ) 


n X-ma 


i-i 


Nidx 


( 10 ) 


L J x 0 


AM em and AMcd are plotted in Figure 2 for each ion as functions of N est . Both 
abundance measures contain fewer ions than the input number, but the regimes of reasonable 
validity for the two measures are different and depend on N cr it- Fortunately, much of the 
discrepancy with the N e or N% volume integrals cancels when we take ratios of the EM or 
CD abundance measures (equivalent to equations (4) and (7), respectively), as shown in 
Figure 3. These plots suggest an approach for the intermediate cases where N est ~ N CT i t . 


First, there must be measurements of N e from (at least) two different density-sensitive 
line pairs: N ei and N e2 . Then, if the column density of the more density-sensitive ion of one 
pair is the numerator and that of the less density-sensitive ion from the other pair is the 
denominator, the column density ratio is greater than the true value and the /^-weighted 
ionic ratio is less than the true value. This follows directly from equations (2)-(7), from 
which we derive 

CD\ _ < Nj, > Net 
CD 2 < N i2 > N ei ' 

Likewise, if the less density sensitive ion of one pair is the numerator and the more density 
sensitive ion from the other pair is the denominator, the column density ratio is less than 
the true value and the iV 2 -weighted ionic ratio is greater than the true value. 
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We note that because N^n for the N III 57 pm line is intermediate between the values 
of N cr i t for the 0 III 52 and 88 pm lines, the estimated N ++ /0 ++ ratio has little uncertainty 
for all values of N e . 

Looking at Fig. 3, we define the three regimes as follows: The low-density regime occurs 
when the estimated density is lower than the lowest N crit for the lines involved (see Fig. 2 
for Ncrit) — here we use the equations for emissivity oc N% since the ionic ratios weighted 
by Nl are close to the model input ratio (unity for the plots). The high-density regime 
occurs where the column density ratio is closer to the model input ratio than the ionic 
ratios weighted by N% — this occurs when the estimated density is higher than the larger 
N C rit of the more density-sensitive line pair. At intermediate densities, both the low-density 
and high-density equations give abundances discrepant from the model input. Considering 
this, in the intermediate regime we average the two cases and quote an uncertainty equal to 
the greater of one half the difference between the two abundance ratios and the statistical 
error. (If one knew for sure that the volumes containing the different ions were co-extensive, 
then an elaborate weighting scheme could be devised, but the overall uncertainties do not 
warrant such schemes.) For the low- and high-density regimes, we will use the equations for 
the ionic ratios weighted by N% or the column density ratios, respectively. Unfortunately, 
it is extremely difficult to get a good estimate of the column density of H + because the 
emissivities of H lines in H II regions and PNe are always proportional to N% and the 
properly weighted density is impossible to determine. Because of this, we are not able to 
measure abundance ratios with respect to hydrogen for all the ions; instead we will obtain 
standard “ionic ratios” with respect to H + only for those ions whose lines have high N^t 
compared to the densities estimated from the observed line pairs. We will calculate both 
ionic ratios and column density ratios for the rest of the heavy element ions observed at FIR 
wavelengths and consider the density regime according to the line’s Ncrit for estimating the 
final abundance ratio. 

A significant difference between this method and previous abundance determinations 
from FIR lines is the increase in the size of the errors that occur when more than one 
density-sensitive line pair is measured in a high-density, compact H II region. The smaller 
statistical errors that one obtains when only one density line pair is measured are an illusion. 
It is clear that measuring multiple density-sensitive line pairs is not just desirable but is, 
in fact, required to obtain reliable abundance measurements from FIR lines (even though 
they may have large assigned errors). We note that a similar situation occurs in optical 
abundance determinations from forbidden lines in that T e must be inferred from various line 
ratios to obtain reliable abundances (e.g., Kennicutt, Bresolin, & Garnett 2003). 
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4. G333.6— 0.2 and W43 

4.1. Analytical Abundance Estimates 

The complete set of derived N e values and abundance ratios is given in Table 2. The 
values of T e used were 6200 K (Caswell & Haynes 1987) and 6500 K (Wink, Wilson, & 
Bieging 1982) for G333.6-0.2 and W43, respectively (mid-and FIR lines are not sensitive 
to uncertainties in T e ). The radio fluxes for G333.6-0.2 were estimated by integrating over 
the model described in the next section at the appropriate distances from the center (Figure 
4) and the radio fluxes for W43 were taken from Lester et al. (1985) (the NRAO VLA Sky 
Survey, Condon et al. 1998, has lower fluxes, probably because the bright clumps in W43 are 
somewhat optically thick at 1.4 GHz, and the “snapshot” observation of Balser et al. 1995 
is missing extended flux). The differences between the ionic abundance results presented 
here and in SCREH for the KAO data are due to different extinction corrections (§2.3), 
different radio fluxes for G333.6-0.2, and different collisional excitation cross sections for 
[S III] (Tayal h Gupta 1999). The differences between the ISO and KAO abundances with 
respect to H + probably indicate the uncertainties in the radio flux estimates. The N ++ /0 ++ 
ratios are the same to within the errors (3a), although some of the difference may be due to 
different beam sizes or pointing errors. 

Estimated abundance ratios for N/S and S/H are given in Table 3, where we assumed 
that the ionization fraction for sulfur <S ++ /S>= 0.9 (a typical value for the low-to-moderate 
excitation and high metallicity models in Rubin 1985, SCREH, and the models computed 
in §5) and that the ionization fraction of <N 3+ /N> is so small that N 3+ can be neglected. 
The N/H ratio was derived by multiplying N/S times S/H. Abundance estimates for the 
other elements are more uncertain because they require the use of models to estimate the 
necessary ICFs. 


4.2. Models of G333.6-0.2 and W43 

We computed H II region models using the program NEBULA (Rubin 1968a; Rubin 
1985; SCREH). We have upgraded NEBULA since the discussion in SCREH in the following 
ways: 

(1) The photoionization and recombination cross sections were extensively updated to 
include results from the Opacity Project (Seaton et al. 1992). Fits to photoionization cross 
sections by Verner et al. (1996) were used for photoionization of C + , C ++ , N°, N + , N ++ , 
0°, 0 + , Ne°, Ne + , S + , S ++ , S 3+ , and Si + . The changes to the coefficients for ionization 
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and recombination of iron are discussed by Rodriguez & Rubin (2004). (Changes to other 
elements were omitted because of the lack of corresponding Opacity Project recombination 
coefficients.) Total recombination coefficients (including dielectronic recombination) were 
taken from the tabulations of Nahar & Pradhan (1997) for C + , C ++ , and N ++ ; Kisielius 
& Storey (2002) for N + ; Nahar (1999) for 0° and 0 + ; Kisielius et al. (1998) for Ne + ; 
Nahar (2000) for S ++ ; and Nahar (1995) for S + and Si + . Charge exchange rate coefficients 
for reactions with H + were taken from the tabulation of Kingdon & Ferland (1996) with 
updates from the Oak Ridge National Laboratory web site “Charge Transfer Database” 3 and 
coefficients for reactions with He + were taken from J. G. Wang et al. (2003, in preparation), 
with the data available from the same Oak Ridge web site. 

(2) The collisional excitation cross sections now include N + and 0 ++ from Lennon & 
Burke (1999), N ++ from Blum & Pradhan (1992), Ne + from Saraph & Tully (1994), Ne ++ 
from McLaughlin & Bell (2000), S ++ from Tayal & Gupta (1999), S 3+ from Saraph & Storey 
(1999) or Tayal (2000) for T e > 10,000 K, Ar + from Pelen & Berrington (1995), and Ar ++ 
from Galavis, Mendoza, & Zeippen (1995, 1998). 

(3) The diffuse radiation field calculation (Rubin 1968a) was modified to include calcu- 
lations at all radial grid positions instead of interpolation between selected positions. Up to 
20 iterations were used if needed to obtain uniform convergence (most positions converged 
within 10 iterations). 

(4) Non-LTE stellar atmospheres with winds from Pauldrach, Hoffmann, & Lennon 
(2001) were used. Smith et al. (2002) have also computed model atmospheres using Paul- 
drach’s code. They note that models computed with these atmospheres give much better 
results for computations of starburst galaxy spectra than the atmospheres from Schaerer & 
de Koter (1997), which were used in the H II region models of Stasinska & Schaerer (1997) 
(we did not use Smith et al.’s models because of the coarseness of their wavelength grid). 
Lumsden et al. (2003) also show that the Pauldrach et al. (2001) atmospheres give better 
results for comparison of measurements of helium lines in Galactic H II regions. 

We used the set of model atmospheres, which were downloaded from Pauldrach’s web 
site 4 , consisting of models for Solar composition, dwarfs (D) and supergiants (S), and T e ff 
= 30, 35, 40, and 45 xlO 3 K. Although G333.6— 0.2 and W43 have metallicities somewhat 
larger than that of the Orion Nebula, the Orion Nebula abundances (e.g., Simpson et al. 
1998) are sufficiently below the “Solar” composition used in the stellar atmosphere models 

3 http : / / ww w-cfadc .phy.ornl. gov / astro/ps/data/ 

4 http : / / w w w . usm . uni-muenchen . de/people/ adi/adi.html 
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that G333.6— 0.2 and W43 are probably closer in metallicity to the “Solar” composition 
models than to anything greatly different (recent work finds a lower Solar O/H ratio, as will 
be discussed later). Stellar atmospheres with higher metallicities have more winds, which 
produce more extreme ultraviolet (EUV) flux, but also more line blanketing (less EUV); 
the difference is probably mainly for frequencies > 4 Ryd, which are not important for H II 
regions. We computed models for W43 using both D and S models from Pauldrach et al. 
(2001) and also the non-LTE, plane-parallel model atmospheres without winds calculated by 
Lanz & Hubeny (2003) 5 . We used only the dwarf star model atmospheres for G333.6— 0.2 
because it is clearly excited by a very young star cluster (Blum et al. 2001), 

One result of the changes to the nebular atomic physics data is that the volume-averaged 
ionization fractions (SCREH; Rubin et al. 1994) for <0 ++ /0> and <Ne ++ /Ne> predicted 
by the models are larger than previous predictions, whereas the fractions for <S ++ /S> and 
<N ++ /N> are approximately the same. The differences are similar to typical uncertainties 
in ionization ratios computed from observational data (< 30%), which would make testing 
difficult unless the data sample was large and of very high quality. However, the choice of 
stellar atmosphere model makes a much larger difference in the ionization fractions and ratios 
thereof; this choice has serious implications for estimates of T e ff and ICFs. We will compare 
various volume-averaged ionization fractions computed for different model H II regions as a 
function of stellar atmosphere model and discuss the need for winds further in §5.1. 


4.2.1. G333.6-0.2 

The initial model is that of a spherically symmetric, core-envelope H II region (Rubin, 
Hollenbach, & Erickson 1993), with density decreasing from high values in the core to low 
values in the outer envelope, as suggested by the estimated densities of Table 2. Our initial 
density distribution for the G333.6-0.2 model computed here was derived by inverting the 
radio brightness temperature distribution, since the radio surface brightness is proportional 
to f Ngdl. The radio brightness temperature distribution of G333.6— 0.2, of necessity, has 
been estimated by combining a variety of different measurements, none of which includes 
the whole H II region at a frequency high enough to insure that the emission is everywhere 
optically thin. The highest spatial resolution radio maps in the literature consist of the 
following maps: (1) a Parkes 64-m telescope map at 8.87 GHz with half-power beamwidth, 
HPBW = 2.5', made by McGee, Newton, & Butler (1979); because the map was made with 
a single dish, it covers the full spatial extent of G333.6— 0.2 to a radius of > 15'; (2) a Fleurs 


5 Available from http://tlusty.gsfc.nasa.gov 
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Synthesis Telescope map at 1.415 GHz with HPBW ~ 50" by Retallack & Goss (1980); only 
the inner 5' could be mapped, and the core is optically thick; (3) an Australia Telescope 
Compact Array map at 3.4 cm (8.8 GHz) with resolution 1.6" by Fujiyoshi (1999); only the 
inner 13" could be mapped, and the peak brightness temperature of ~ 4400 K indicates that 
the central 2" — 3" is not completely optically thin (T e ~ 6200 K, Cas-well &; Haynes 1987). 
However, this last map is extremely important because it shows just how compact and bright 
the central core of G333.6— 0.2 really is. A compact bright core is also required by the NIR 
maps and the density-sensitive FIR line ratios. 

Still lacking is a measurement of the transition region between the inner core and the 
envelope. There is, however, a map of the warm dust emission made by the Midcourse Space 
Experiment (MSX; Price et al. 2001) at 21 pm with 19" resolution. 6 Although there is no 
certainty that warm dust has the same spatial structure as ionized gas (and it certainly does 
not have constant temperature), in the absence of an alternative -we used it to interpolate 
between the core and envelope radio maps. 

The total radio flux was taken from the measurements at the highest observed fre- 
quencies because G333.6— 0.2 shows evidence of being at least partially optically thick at 
frequencies < 10 GHz. We used the integrated flux of 93.8 Jy at 14.7 GHz measured by 
McGee & Newton (1981). The final estimated radio brightness distribution was renormal- 
ized to produce this total integrated radio flux. The estimated radio brightness temperature 
distribution, scaled oc u~ 21 as though it were optically thin at 5 GHz, and derived electron 
densities are shown in Figure 4. The total radio flux requires a stellar ionizing photon lumi- 
nosity of 10 5 ° photons s -1 . There may also be, of course, additional EUV photons that do 
not contribute to the H II region photoionization because they are absorbed by dust. Table 4 
lists the relevant IR and radio observations, corrected for extinction, including those mid-IR 
observations from the literature which include most of the core. 

In addition to matching the radio brightness temperature distribution, the observations 
of G333.6— 0.2 put a number of requirements on the parameters of any model. There must 
be a core of very high density to produce the bright radio and mid-IR core and the high 
[O III] 52/88 pm and [S III] 19/33 pm line ratios seen in the center position. The material 
at 2' — 4' must be of low density to produce the 52/88 pm and 19/33 pm line ratios of order 
unity (N e of order 150 — 400 cm -3 ). The material further out must be of even lower density 
(7V e ~ 40 cm -3 ) to make a large contribution to the 122/205 pm line ratio ~ 2.5. The 
stellar T eff has to be high enough to produce 0 ++ and N ++ at all positions (out to > 3') but 
low enough so that there is ample N + , Ar + , and not much N ++ , 0 ++ , or S 3+ even in the 


6 Available at http://irsa.ipac.caltech.edu 
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center position. In addition, the He + /H + radio recombination line (RRL) ratio is only 0.045 
(McGee & Newton 1981), indicative of a relatively low T e n (< 35,000 K). 

We computed over a hundred models of G333.6— 0.2 with a variety of combinations of 
stellar atmosphere models and density filling factors, all models reproducing the observed 
radio flux. The abundance ratios used in the models were estimated from the fluxes in 
Table 4 (see also Table 3) and range from 1.25 to 3.4 times the Orion Nebula abundances 
(Simpson et al. 1998). Line fluxes were" computed from the models by integrating over 
flat-topped circular beam profiles for the IR lines or circular Gaussian beam profiles for the 
radio recombination lines located at the projected distances of the observed positions from 
the H II region center. 

Our best example of a spherically-symmetric model is the G333.6— 0.2 Core model in 
Table 4. This particular model is shown because it predicts line ratios that agree reasonably 
well with the density indicators ([S III] 19/33 pm, [O III] 52/88 pm, and [N II] 122/205 
pm) and the ionization indicator measured in the largest beam (He + /H + RRLs measured 
by McGee & Newton 1981). The model uses a linear combination of 0.9375 times the D-30 
(30,000 K Dwarf) plus 0.0625 times the D-35 stellar atmosphere models of Pauldrach et 
al. (2001) and has a volume filling factor of 0.5. The contribution from the D-35 model, 
although seemingly small, is necessary to have enough high energy photons to ionize He and 
N to the observed ionization. The given model combination only represents the shape of the 
ionizing spectrum and not the actual numbers or spectral types of the ionizing stars. 

However, no model that we computed is a good fit — the data seem to require stellar 
ionizing fluxes with both higher and lower average T e ff- The large fluxes of low-excitation 
ions, such as Ne + and Ar + in the center of G333.6-0.2 seem to require stars with low T eff 
(< 32,000 I< was suggested by Rank et al. 1978). However, the central density and optical 
depth are so high that when the exciting star is this cool, all the high energy photons that 
are capable of ionizing 0 + to 0 ++ are absorbed at small radii (0.87 pc, corresponding to 
60") and the models have essentially no 0 ++ at the radii corresponding to our three North 
positions. Yet we observe 0 ++ at all four positions. It would seem that an exciting star with 
a somewhat higher T eff (> 35, 000 K) is required or a combination of stars including a much 
hotter star; however, for such models the > 35 eV photons required to ionize 0 + to 0 ++ still 
are confined in the central core but the total [O III] line fluxes are much larger than observed 
unless the O/H, O/N, O/Ne, etc. ratios are unreasonably low. On the other hand, in H II 
regions ionized by stars with relatively low T e g, the N ++ Stromgren radius, Rs, is much 
larger than the 0 ++ Rs, leading to observed N ++ /0 ++ ratios that can be much larger than 
the actual N/O ratio (e.g., Rubin et al. 1994; SCREH; Stasinska & Schaerer 1997). The 
result of the differences in the N ++ and 0 ++ Rs is that the observed N ++ /0 ++ ratio should 
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increase with projected distance from the center. This predicted increase is not observed in 
our N ++ and 0 ++ data taken at the North positions. Likewise, the ISO data in a 80" beam 
taken at the center position do not produce a larger N ++ /0' t ' + ratio than does our KAO 
data in a 45" beam. Again, the observed behavior of this line ratio can be reproduced only 
by models with higher T e g stars. In fact, all the models that we computed predict a much 
larger [N III] 57/[0 III] 52 /im ratio than is observed at any position. Agreement cannot be 
improved by increasing the O/H ratio (the models already predict more [O III] 52 /zm flux 
than is observed) or decreasing the N/H ratio (the model [N II] lines fluxes agree with the 
observations). A lower T e ff seems to be required, but already there is too little 0 ++ in the 
North positions. We conclude that a spherical model, even of varying density, cannot satisfy 
the observations of both the very bright core and the doubly ionized lines in the low density 
region at 2' — 4' from the center. 

A possible solution might be some sort of non-spherical model with large clumps, or a 
blister model (face-on from its symmetrical appearance, Hyland et al. 1980). The NEBULA 
code for a blister model (Rubin et al. 1991) assumes that the density distribution is that of 
the surface of a molecular cloud with some possible gradient of density going into the cloud. 
The stellar photons ionize a hole on the edge of the cloud. However, that axisymmetric model 
does not use a density function that allows for a particularly high density immediately next 
to the stars, as is required for a model of G333.6— 0.2. Without a modification to the 
axisymmetric 2-D code, this code is not suitable. A treatment that can handle a two- or 
three-dimensional density distribution appears to be required. 

We next considered a quasi-blister model or two-component model, such as was used 
by Simpson et al. (1986) and Colgan et al. (1991) to model the Orion Nebula and K3-50, 
respectively. This would represent either a blister H II region with one open side or an H II 
region with very large clumps where ionizing photons escape through holes or paths between 
the clumps. A linear combination was made of two spherical models, one with the high 
density core and one with an almost empty center to provide the doubly ionized ions seen 
at 2' - 4'. The model with the empty center had a very large Rs and therefore had to be 
arbitrarily truncated in order that it be density bounded at about the same radius as the 
Rs for the Core model. An example of such a shell model, labeled “Model for the North 
Positions” , is also given in Table 4; the model used the same stellar atmosphere as the Core 
model but the filling factor was unity. Even here there are not enough photons to ionize 0 ++ 
at the N2 and N3 positions and too much N ++ at the Nl position. Possibly three or even four 
spherical models with different interior shell diameters and different scaling fractions for the 
volume of the model in the line of sight are needed to reproduce the observations at the four 
observed positions. Because of the arbitrariness of the procedure and the non-uniqueness of 
the result, we did not calculate any of these multi-shell models. 
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The abundances have been only partially adjusted to produce the observed fluxes — the 
immediate goal of the models was to determine T e g, the filling factor, and other geometry 
factors. Further adjustment to the abundances is not meaningful because the geometry 
still is not well determined, although the agreement of the He + /H + RRL ratios with the 
observations indicates that T e ff is not especially discrepant (since most helium is primordial, 
the interstellar He/H ratio is not increased significantly as a result of Galactic chemical 
evolution). The uncertainty in the geometry factors makes it difficult to estimate a reliable 
abundance for oxygen. However, if we consider just the Center and N1 positions, we find 
that the observed fluxes would be matched by the model fluxes if 0/H~ 3 x 10 -4 , a value 
somewhat lower than that observed in the Orion Nebula. This O/H abundance ratio is not 
reliable, but suggests that O/H is not especially high in G333.6— 0.2. As expected, since we 
observed the dominant ionization states of nitrogen and sulfur (N + and S ++ ) in this source, 
our final estimated abundances of these elements are those input to the models (Table 3). 
On the other hand, it is difficult to think of a geometry that would produce the large [Ar II] 
6.98 /rm line flux observed in the G333.6— 0.2 core, given the number of high energy photons 
that are needed to produce the observed [Ne III] and [O III] lines. This discrepancy may 
indicate that the atomic data for argon need revision (see also Morisset et al. 2003). 

Kim &; Koo (2001) showed that many, if not most, ultracompact H II regions are not 
ionization bounded at the radius given by the ultracompact component, but have extended 
envelopes ionized by the same stars. They suggested that the ultracompact components are 
blister H II regions and the extended envelopes are the result of champagne flows combined 
with hierarchical clumping. Although G333.6— 0.2 is larger than the typical ultracompact 
H II region because of its multiple stars and high ionizing luminosity, it does have the 
high central density and optical depth of an ultracompact H II region. We conclude that 
G333.6— 0.2 may be an example of a large ultracompact H II region with an extended enve- 
lope. 


4.2.2. W43 

Because the constant density model by SCREH for W43 is not a good match to the 
densities inferred from the observations (see Table 2), we calculated new models for it. 
The appearance of W43 is that of a large “T” or incomplete “D” shape projected on the 
sky (Lester et al. 1985; Liszt 1995; Two Micron All Sky Survey image 7 ) with a central 
star cluster near but not on the vertical element of the letter. Because our three observed 


7 2MASS, http://irsa.ipac.caltech.edu 
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positions, including the one closest to the cluster, all have similar densities and excitation 
(SCREH; Table 2), we interpret the geometry as a cluster of ionizing stars surrounded by 
a cavity and then a shell of ionized gas, such that all three positions are at essentially the 
same physical distance from the stars. Thus we model the H II region as a shell of ionized 
gas with interior radius = 1.2 pc (0.75' at a distance of 5.6 kpc (Reifenstein et. al. 1970; Liszt 
1995), scaled to a Galactic Center distance of 8 kpc). The model N e drops smoothly from 
770 cm -3 at 1.2 pc to 420 cm -3 at Rs ~ 2.9 pc, matching the average of the total observed 
flux of 77.5 Jy at 2 cm (Schraml & Mezger 1969) and 97.4 Jy at 6 cm (Reifenstein et al. 
1970) for an ionizing luminosity of 2.9 x 10°° photons s -1 . The densities and ionic abundance 
ratios of the three positions were averaged (weighted) for comparison with the models; these 
comparisons are given in Table 5. Three models were computed with the only difference 
being the input stellar atmosphere models, as described in the Table. The reason the S 
(supergiant) model atmospheres were also used is the presence in W43 of a Wolf-Rayet star 
and two O giants/supergiants (Blum et al. 1999). The combinations of stellar atmosphere 
models were adjusted until the H II region models produced the observed average radio 
He + /H + ratio (He/H = 0.10 for the nebular models) and the observed N ++ /N + ratio. 

To correct the observed abundance ratios for unseen ionization states, we used the ion- 
ization fractions calculated from the model as shown in the Table 5 comments and foot- 
notes. The biggest difference between the models is for the <0 ++ /0>/<S ++ /S> and 
<N ++ /N>/<0' f+ /0> ratios, no doubt because the Supergiant model atmospheres have 
less flux from 35 to 41 eV than the Dwarf model atmospheres. Pauldrach et al.’s (2001) 
Supergiant models have much higher mass loss rates than their Dwarf models, so the differ- 
ence must be caused by more line blanketing at these energies by the gas in the wind. On 
the other hand, the static, plane-parallel, non-LTE atmospheres of Lanz &: Hubeny (2003), 
which have no winds, have much lower fluxes from ~ 40 to ~ 50 eV. 

All three W43 models adequately reproduce the T e ff-indicating indices: the He + /H + and 
N ++ /N + ratios. The range of implied abundances in Table 5, though, indicates the large 
uncertainties resulting from determination of abundances through model-produced ICFs, 
particularly when different stellar atmosphere models are used to compute the ICFs. We 
prefer the values of N/H derived from comparisons with sulfur, rather than those derived 
from oxygen. The estimated N/O ratio is 0.32 to 0.42 from ICFs computed with models D 
and S (not very different from the 0.34 estimated by SCREH), but is as low as 0.16, similar 
to that of the Orion Nebula (Rubin et al. 1991), from the ICFs computed with model T. 
The reason for the wide range is not the abundances of nitrogen or sulfur, which are reliable 
because they do not require large ICFs. The problem is the total abundance of ox 3 >-gen: the 
model ICFs produce O/H ranging from 2.8 to 8.3 xlO -4 , whereas the predicted O/H ratio is 
7.2 x 10 -4 for a Galactic abundance gradient of d(log O/H )/dR = —0.06 dex kpc -1 (Henry 
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& Worthy 1999), given the Orion Nebula’s location at Rq ~ 8.4 kpc and O/H = 4.0 to 
4.4 x 10 -4 (Rubin et al. 1991; Esteban et al. 1998). In fact, one expects that the oxygen 
abundance for W43 should be much higher than that of the Orion Nebula, given W43’s low 
T e = 5410 to 6500 I< (Subrahmanyan &; Goss 1996; Wink, Wilson, & Bieging 1983; and 
references in both papers). For this source, as for G333.6-0.2, we conclude that we still do 
not have a good estimate for the abundance of oxygen, the most abundant heavy element 
and the most important contributor to cooling in H II regions. 


5. Discussion and Conclusions 
5.1. Stellar Atmosphere Models 

As was discussed above, many of the ionization correction factors are critically dependent 
on the choice of stellar atmosphere model. This can be demonstrated either by plotting 
ratios of the numbers of photons able to ionize various ions as a function of T e ff or by 
plotting ionization fractions or ICFs calculated from models; here we do the latter because 
ICFs are more readily applicable to abundance computation. This is illustrated in Figure 
5, where we plot the ICF ratios <N ++ /N>/<0 ++ /0>, <0 ++ /0>/<S ++ /S>, and 0 ++ /0 
versus the abundance-independent excitation ratio <N ++ /N>/<N + /N>. The models from 
Rubin (1985) that are plotted in Fig. 5 all have total nucleon density = 10 3 cm -3 , ionizing 
luminosities of 10 49 to 10 5 ° photons s _1 , and Orion Nebula metallicities; the models from 
SCREH had similar densities and ionizing luminosities but metallicities ranging from 1 to 
2 times Orion nebula metallicities. The models from both Rubin (1985) and SCREH were 
computed with the versions of NEBULA available at the time and both used LTE stellar 
atmosphere models from Kurucz (1979). In addition to the models described in §4, constant 
nucleon density (= 10 3 cm -3 ) spherical models with ionizing luminosities of 10 49 photons s' 1 
were computed and plotted. These recent NEBULA models used stellar atmosphere models 
from either Pauldrach et al. (2001), Lanz Sz Hubeny (2003), or contemporary 8 LTE stellar 
atmosphere models described by Kurucz (1992). The stellar atmosphere models had “Solar” 
abundances and the nebular models had 1.5 times Orion Nebula abundances (Simpson et al. 
1998) such that O/H = 6.0 x 10 -4 . Because they are widely used, we also plot ICFs from 
the models of Stasiiiska & Schaerer (1997). For these models both the stellar and nebular 
abundances were “Solar”, with 0/H= 8.51 x 10 -4 (currently it is believed that the Solar 
O/H ratio is lower: 5.45 x 10 -4 ; Holweger 2001). 


8 http: //kurucz. harvard.edu/ 
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Unless an H II region can be observed at both optical and infrared wavelengt hs, observers 
cannot measure the abundances of N, 0, Ne, S, and Ar without making corrections for 
unseen ionization states for some elements. For optically obscured H II regions like W43 
and G333.6— 0.2, the most significant missing ion is 0 + , since oxygen is the most abundant 
heavy element and since these low excitation H II regions are believed to have most of their 
oxygen as 0 + . Rubin et al. (1994) and SCREH assumed that they knew the O/S ratio and 
used the observed 0 ++ /S ++ ratio to estimate the ICF for their measured N ++ /0 ++ ratios, 
since they did not have any other line pair that could be used to indicate excitation. In this 
paper we use our observed N ++ /N + ratios to estimate T^\ however, we see from Fig. 5a 
how dependent the N ++ /0 ++ ICF is on the stellar atmosphere model. 

Martm-Hernandez et al. (2002a) took advantage of their ISO SWS measurements of the 
Ne ++ /Ne + ratios to estimate ICFs without the need for assuming a priori some abundance 
ratio; unfortunately, the large beam size difference between the SWS and LWS instruments 
means that their observed Ne ++ /Ne + ratios may not always be appropriate for estimating 
the excitation of 0 ++ and N ++ . Moreover, Ne ++ requires ionizing photons >41 eV, and it 
is at these high energies that the stellar atmosphere models are especially divergent, as we 
now discuss. 

An important test of the validity of stellar atmosphere models of hot stars is whether 
H II region models produced with these atmospheres predict line fluxes that agree with 
observations. Except for a few rare instances, this is the only way that the EUV fluxes of 
the models can be tested, since interstellar absorption by hydrogen predominantly prevents 
the EUV fluxes from hot 0 stars from reaching the Earth. Here we discuss the stellar fluxes 
that can doubly ionize neon ( hv >. 41 eV). This has been called the [Ne III] problem (e.g., 
Sellmaier et al. 1997) because it has been observed that the Ne ++ /0 ++ ratio is relatively 
constant over a large range of H II region excitation (SCREH; Stasinska & Schaerer 1997 
and references therein; Kennicutt et al. 2003), contrary to predictions of models that used 
LTE stellar atmospheres (e.g., SCREH). Proposed solutions have involved using non-LTE 
stellar atmospheres, usually with winds (Rubin, Kunze, k Yamamoto 1995; Sellmaier et al. 
1996; Stasinska & Schaerer 1997). Figure 6 shows the observations of SCREH along with 
the predictions of the nebular models discussed in this paper. Models calculated with black 
bodies for the stellar atmosphere spectrum are also plotted since Morisset et al. (2003) found 
that such models give a reasonable fit to the mid-IR ISO observations. 

Surprisingly, in Fig. 6 only the models of Stasinska k, Schaerer (1997) reproduce the 
observations over the total range of observed 0 ++ /S ++ , although if Ne/O is as large as 0.25 
in all H II regions as it is in the Orion Nebula (Simpson et al. 1998; Table 3), these H II 
region models produce too high a Ne ++ /0 ++ ratio because of the very large EUV fluxes in 
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Schaerer & de Koter’s (1997) stellar atmosphere models. Moreover, Schaerer & de Koter’s 
models have too much EUV flux for reproducing other observations (e.g., Smith et al. 2002). 
Earlier, Sellmaier et al. (1997) had demonstrated that the [Ne III] problem was supposedly 
solved when they obtained a good fit to the data by using non-LTE atmospheres with winds 
computed with Pauldrach’s code as it then existed. However, our H II region models with 
non-LTE stellar atmospheres with winds from Pauldrach et al. (2001) taken from both A. 
Pauldrach’s web site and from the stellar atmosphere models produced by Smith et al. (2002) 
using Pauldrach’s WM-basic code predict much lower Ne ++ /0 ++ than observed for Dwarf 
atmospheres with T e ff < 40, 000 K and for Supergiant atmospheres with T e ff < 35, 000 K. 
Our H II region model produced with a 35,000 K supergiant atmosphere from F. Sellmaier 
(private communication in 1995) lies above the Pauldrach et al. Supergiant line but not as 
high as the model with the 35,000 K supergiant atmosphere in Sellmaier et al.’s (1996) Figure 
2, even though we are now assuming Ne/O = 0.25 (Simpson et al. 1998) instead of 0.2025. 
Part of the difference of the current nebular model predictions with previous models is the 
higher 0 ++ ionization relative to Ne ++ and S ++ obtained with the Opacity Project cross 
sections described in Section 4. The model computed with Sellmaier’s (1995) atmosphere has 
higher ionization than the model computed with Pauldrach et al.’s (2001) S-35 atmosphere 
because it used Orion Nebula abundances (as in Sellmaier et al. 1996) instead of the higher 
metallicity (approximately 2 times the Orion Nebula abundances of Table 3) used for the rest 
of the models plotted here. The result, though, is a return of the mismatch between the stellar 
atmosphere models and the EUV fluxes of real stars, with possible additional discrepancies 
for the ionization, recombination, and collisional excitation cross sections for oxygen, neon, 
and sulfur. Moreover, we see that all the models computed without winds do a poorer job 
of reproducing the observations than both the models computed with non-LTE atmospheres 
with winds and the models computed with simple black bodies for stellar atmospheres (for 
this reason we give low weight to the T model for W43 in Table 5). This probably means 
that the codes producing the stellar atmosphere models used for computation of H II region 
models should include the effects of stellar winds. 


5.2. Abundances and Abundance Gradients 

At the low excitations seen in inner Galaxy H II regions, the predicted flux of the 0 ++ 
lines is extremely sensitive to the stellar atmosphere, but if the models are accurate, the 
O/H ratio of both G333.6— 0.2 and W43 could be as low as 3 — 5 x 10 -4 , substantially 
lower than would be predicted from the observed Orion Nebula abundance and an O/H 
abundance gradient of —0.06 dex kpc -1 (Henry & Worthey 1999). This low O/H ratio 
is hard to understand because of the low T e measured in inner Galaxy H II regions from 


-23 - 


radio recombination lines. There would need to be some source of systematic error for 
the O/H ratio to be higher than these estimates: i.e., what would be required is that the 
excitation of oxygen in the two H II regions is lower than estimated in this paper from 
He + /H + recombination lines and N + and N ++ forbidden lines. Possible reasons that the 
G333.6— 0.2 and W43 excitation might be overestimated are as follows: (1) The abundance 
of N + is underestimated because the [N II] line fluxes are undermeasured owing to telluric 
absorption at 122 /xm or diffraction at 205 /xm. This is a particular problem for W43, where 
the two measurements of the [N II] 122 /xm line in different apertures can be used to produce 
very different electron densities and N ++ /N + abundance ratios (1.5 to 3.0 in Table 5). (2) 
The estimated N ++ /N + ratio is too high because some of the estimated N + is missing owing 
to inadequate correction jor density variations. (3) The estimated N ++ /N + ratio is too 
high because of incorrect collision al excitation cross sections for the pertinent energy levels. 
(4) The models predict too much 0 ++ /0 versus N ++ /N + . The spread of 0 ++ /0 ratios 
in Fig. 5c shows that any excessive 0 ++ /0 is due to model atmospheres and not due to 
the ~ 30% higher oxygen ionization resulting from the Opacity Project cross sections. (5) 
The ionization of helium might be overestimated if the He/H abundance ratio is actually 
greater than the assumed 0.10. We conclude that determinations of the O/H abundance 
ratio in low excitation H II regions have large uncertainties at this time, with the largest 
contributor to the uncertainty being the choice of stellar atmosphere models, with some 
possible contribution from the atomic data used in the nebular models. 

We have demonstrated that' the abundances of nitrogen and sulfur can be obtained from 
FIR observations and the ratio of N/S can be obtained when both ionization states of N + and 
N ++ are measured. (Simpson et al. 1998, Martfn-Hernandez et al. 2002a, and Giveon et al. 
2002 also show that the abundances of Ne/H and Ar/H, as well as S/H, can be determined 
from mid-IR observations with little bias owing to excitation or abundance.) Compared to 
SCREH, the S/H abundance ratio is higher in G333.6— 0.2, owing to the new data at the 
additional North positions, but lower in W43 because of the higher collisional excitation 
cross sections for sulfur. The N/H ratio is about 50% larger in G333.6— 0.2, owing to the 
new data at the North positions, but a factor of 2 smaller in W43, owing to the smaller ICF 
for N + . However, since the estimated O/H ratio is also much smaller in W43, now that we 
are estimating the ionization from N ++ /N + and not from 0 ++ /S ++ with the assumption 
of a constant O/S ratio as SCREH did, the estimated N/O ratios agree with SCREH for 
both G333.6— 0.2 and W43. Thus we conclude that with this new analysis, the S/H and 
N/O ratios are consistent with the abundance gradients measured by SCREH, Rudolph et 
al. (1997), and Afflerbach et al. (1997), and Martfn-Hernandez et al. (2002a) for N/O. The 
N/H ratio for G333.6— 0.2 is also consistent with the gradient of —0.10 dex kpc -1 estimated 
by SCREH and Afflerbach et al. (1997) but the N/H ratio inferred for W43 is lower than 
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what would be expected from the previously published abundance gradients. 

On the other hand, the estimated O/H ratio is consistent with no gradient at all, 
working inward from the Orion Nebula. This is in contradiction to the other abundance 
studies of H II regions and PNe in the Galaxy, particularly in regions other than the inner 
Galaxy where optical studies allow measurement of 0 + (see the reviews of Henry & Worthy 
1999 and Rolleston et al. 2000 and references therein). Smartt et al. (2001) measured 
abundances in four B stars in the inner Galaxy (Rq ~ 2.5 to 5 kpc) by optical absorption 
line spectroscopy. They found that these four stars have oxygen abundances similar to 
stars in the Solar neighborhood, even though they have enhanced abundances of nitrogen, 
magnesium, aluminum, silicon, and sulfur corresponding to their inner Galaxy locations 
and expected abundance gradients of —0.07 dex kpc -1 . They have no explanation for this 
anomalous oxygen abundance, nor do we. We conclude that more measurements of the 
oxygen abundance in inner Galaxy stars and H II regions are needed to confirm this apparent 
lack of increase in the O/H ratio, because we know that oxygen is one of the main coolants 
in H II regions and the low observed radio T e measurements require additional cooling from 
heavy elements compared to Solar-neighborhood H II regions. 


5.3. Conclusions Regarding Systematic Errors 

There are several possible sources of systematic error that must be considered when 
measuring abundances in obscured, low-excitation H II regions, that is, all the H II regions 
in the inner Galaxy. The impact of these systematic errors is that there is a very large 
uncertainty in the total metallicity and important abundance ratios like N/O. These sources 
of systematic error are the following: 

(1) The most abundant heavy element ion in low-excitation H II regions, 0 + , has no 
bright infrared lines. The consequence is that the abundance of this critically important ele- 
ment must be estimated from the abundance of an ionization state with fractional abundance 
< 0.5. 

(2) Nitrogen, the secondary element which one needs for studies of galaxy chemical 
evolution, has two important ionization states: N + as well as N ++ . Fortunately, both species 
have measurable FIR lines, but both must be measured to determine the excitation. 

(3) Moreover, both FIR N + lines have greatly different values of from those of 
S ++ or N ++ . Thus one needs to be aware of and compensate for the ever-present density 
variations. 
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(4) The average density could be quite different in the various ionization zones, as shown 
by the different electron densities derived from A^-diagnostic line pairs: low ionization (N + ), 
intermediate ionization (S ++ ), and high ionization (0 ++ ). 

(5) ICFs estimated from models have major uncertainties owing to uncertainties in the 
stellar atmosphere models and possibly the atomic data used in the H II region models. 

It is clear that the abundances of primary (oxygen, neon, or sulfur) and secondary 
(nitrogen) elements will need to be measured in many more inner Galaxy H II regions before 
we understand the chemical evolution of the Milky Way. In the near future, measurements 
of N ++ , N + , 0 ++ , Ne ++ , and S ++ will be possible from the Stratospheric Observatory 
for Infrared Astronomy, SOFIA. We recommend that future observers measure excitation- 
sensitive lines, such as Ne + and Ne ++ , as well as both N + and N ++ to determine the H II 
region ionization state needed to estimate the ICF for 0 + . It is important that a sufficient 
number of lines be measured in each H II region (and extended H II regions mapped) so 
that detailed models of the H II region can be made. Discrepancies between the predicted 
line fluxes from the models and the observations can then be used to indicate EUV energy 
regions where the stellar atmosphere models may need revision. 
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Fig. 1. — (a) Ratios of volume emissivities for line pairs of the electron-density-sensitive 
lines [S III] 19/33 ^ra, [0 III] 52/88 pm, and [N II] 122/205 pm are plotted as a function 
of electron density, N e , for electron temperature, T e , of 7000 K. (6) Line intensities were 
calculated for slab models with N e — N^e~ x by integrating the volume emissivity of each 
line from x 0 to x max . N 0 ranges from 1 to 10 6 cm -3 and N e = 0.01 cm -3 at x max (see text). 
Values of N e were estimated from the abscissa of Fig. la by assuming that the ratios of the 
line intensities could be used for the ordinate. These estimated values of N e divided by Nq 
are plotted versus N 0 . The root-mean-squared (RMS) density divided by N 0 is also plotted 
(dotted line). 

Fig. 2. — Line intensities were calculated for slab models with N e = N 0 e~ x by integrating 
the volume emissivity of each line from x-o to x max \ corresponding densities were estimated 
from the density-sensitive line ratios in Fig. la. The number of ions estimated from the 
line fluxes and estimated volume emissivities are plotted for lines of sight with varying 
amounts of high-density gas (see text); the abscissa is the estimated N e for the line of sight. 
These abundance measures, AMem and AMqd > are defined as the calculated line intensity 
divided by the normalized volume emissivity or the emissivity per ion, respectively. They are 
normalized by dividing by the emission measure or by the column density of ions, respectively 
(see text). Because preferred abundance measures for each ion depend on the density, critical 
densities, N crit , are also tabulated for T e — 7000 K. (a) S ++ measured from [S III] 19 and 
33 pm lines. ( b ) 0 ++ and N ++ measured from [O III] 52, [O III] 88, and [N III] 57 pm lines. 
Note the near coincidence of the N ++ and 0 ++ abundance measures, (c) N + measured from 
[N II] 122 and 205 pm lines. 

Fig. 3. — Ratios of the abundance measures for the ions shown in Fig. 2 are plotted as a 
function of estimated N e for lines of sight with varying amounts of high-density gas. The 
heavy line is AMem and the dotted line is AMqd (the light line is unity, for reference), 
(a) N ++ /0 ++ . The lines for AMem and AMqd coincide. (6) 0 ++ /S ++ . (c) N + /S ++ . (d) 
N++/S++. ( e ) N ++ /N + . 

Fig. 4. — (a) Azimuthally-averaged radio brightness temperatures for G333.6— 0.2 are plotted 
as a function of angular distance on the sky from the central peak. The radio brightness 
temperatures were all scaled to the same frequency, 5 GHz, using the formulation of Mezger 
& Henderson (1967) and ignoring optical depth effects. ( b ) N e , derived by inverting the radio 
brightness temperature distribution shown in Fig. 4a, is plotted as a function of distance, 
r, from the central star cluster. A spherical H II region at a distance of 3 kpc was assumed. 
The overall form of this density function is N e oc r -12 , as shown by the dashed line, which 
was fit by least squares. 
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Fig. 5. — Ratios of ionization fractions are plotted for spherical H II region models. The 
models with LTE stellar atmosphere models are those from Rubin (1985) and SCREH, 
who used model atmospheres from Kurucz (1979), and the new NEBULA models with 
contemporary atmosphere models described by Kurucz (1992) (r eff = 33, 35, 37, 40, and 
45 x 10 3 K). The H II region models with non-LTE stellar atmosphere models are those 
from Stasinska & Schaerer (1997), who used atmosphere models from Schaerer & de Koter 
(1997), and the NEBULA models calculated for this paper, where we used non-LTE at- 
mosphere models from Pauldrach et al. (2001) and Lanz & Hubeny (2003) (T e ff = 30, 
35, 40, and 45 x 10 3 K). For the H II region models calculated with Pauldrach et al. at- 
mospheres, the solid line connects models with Dwarf atmospheres and the dashed line 
connects models with Supergiant atmospheres. For the H II region models calculated with 
Lanz & Hubeny atmospheres, the solid line connects models with atmospheres with log 
g = 4.0 and the dotted and dashed lines connect models with atmospheres with log g = 3.0 
to 3.5 and with Lyman continuum luminosities of 10 49 and 10 50 photons s _1 , respectively. 
The notation in the legend applies to all parts of the figure, (a) <N ++ /N>/<0 ++ /0> 
versus <N ++ /N>/<N + /N>. ( b ) <0+ + /0>/<S ++ /S> versus <N ++ /N>/<N + /N>. (c) 
<0 ++ /0> versus <N ++ /N>/<N+/N>. 

Fig. 6. — Ratios of Ne ++ /0 ++ inferred from the observations of Simpson et al. (1995) are 
plotted against their inferred 0 ++ /S ++ ratios. Ratios from the nebular models discussed 
in this paper are also plotted. The notations of Pauldrach et al. (2001), Sellmaier (1995), 
Lanz & Hubeny (2003), and Kurucz (1992) in the legend refer to the stellar atmospheres 
used in our NEBULA models (same T e s as the models in Fig. 5; Sellmaier’s atmosphere is 
a Supergiant with T e ff = 35,000 K). The plotted Ne ++ /0 ++ ratio for the models is equal 
to the calculated <Ne ++ /Ne>/<0 ++ /0> ratio multiplied by an assumed abundance ratio 
of Ne/O = 0.25 and the plotted 0 ++ /S ++ ratio for the models is equal to the calculated 
<0 ++ /0>/<S ++ /S> ratio multiplied by an assumed abundance ratio of O/S = 50 (Orion 
Nebula abundance ratios from Simpson et al. 1998). For the H II models calculated with 
Pauldrach et al. atmospheres, the solid line connects models with Dwarf atmospheres and 
the dashed line connects models with Supergiant atmospheres. For the H II region models 
calculated with Lanz & Hubeny atmospheres, the solid line connects models with atmo- 
spheres with log g = 4.0 and the dotted and dashed lines connect models with atmospheres 
with log g = 3.0 to 3.5 and with Lyman continuum luminosities of 10 49 and 10°° photons s -1 , 
respectively. Models computed with atmospheres represented by black bodies from 20,000 K 
to 40,000 K (every 5000 K) are also plotted (open circles connected by a dotted line). 
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Table 1. Observed Line and Continuum Fluxes. 


Line 




Line Flux® (10“ 18 W cm" 2 ) 
Continuum Flux® (Jy) 




G333.6-0.2 


W43 (G30.8-0.0) 

G333.6-ISO b 

G333.6-N1 C 

G333.6-N2 C G333.6-N3 C 

W43N C 

W43C C 

W43S C 

[S III] 18.7 ftm 


160.1*16.2 

36.3*6.4 


SCEEH 

SCREH 

SCREH 



1370*275 

290*280 


SCREH 

SCREH 

SCREH 

[S III] 33.5 /im 


189.1*3.6 

65.7*4.7 

38.9*3.1 

SCREH 

SCREH 

SCREH 



7040*580 

740*185 

730*180 

SCREH 

SCREH 

SCREH 

[0 III] 51.8 /im 

178.2*8.5 

84.3*2.6 

13.1*1.6 

6. 8*1. 2 

SCREH 

SCREH 

SCREH 


54400 

9910*320 

1130*360 

410*270 

SCREH 

SCREH 

SCREH 

[N III] 57.3 fJtn 

46.3 *4.9 

51.1*2.8 

11.9*1.4 


SCREH 

SCREH 

SCREH 


57630 

10945*230 

1310±200 


SCREH 

SCREH 

SCREH 

[0 III] 88.4 /im 

39.9*3.9 

37.7*0.8 

14.5*0.9 

7. 5*0. 7 

SCREH 

SCREH 

SCREH 


55100 

9030*175 

1060*270 

630*90 

SCREH 

SCREH 

SCREH 

[N II] 121.9 /im 


7.4±1.4 

2.6±0.8 


14.4*0.8 

5. 5*0. 5 

7. 9*0.8 



6790*120 

575*230 


5270*420 

1450*260 

4690*380 

[N II] 121.9 /im d 





4.9*0.5 d 

2.5±0.4 d 

2.6±0.6 d 






2530*21 0 d 

800±200 d 

1600*300 d 

[N II] 205.2 /im 


20*0.2 

1. 5*0.2 

1.7*0. 2 

1.52*0.26 

0.46*0.14 

0.56*0.13 



2000*145 

230*120 

115*100 

1300*70 

490*50 

1440*40 


a Flux and continuum errors are 1 a and are statistical only. 
b All ISO LWS data had entrance aperture 80". 
c All CGS data had entrance aperture 60" except as marked. 
d Entrance aperture 45". 
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Table 2. Electron Densities and Abundance Ratios 





G333.6-0.2 




W43 



G333 6-ISO 

G333.6-C a 

G333.6-N1 

G333.6-N2 

G333.6-N3 

W43N 

W43C 

W43S 

7*9.6 

1.5 

1.5 

1.5 

1.5 

1.5 

2.74 

2.74 

3.63 

S'sGHz*’ 

40.75 

26.65 

5.35 

1.10 

0.35 

5.91 

2.25 

3.11 

Electron Densities (cm 3 ) 








Ne (S III) 


2400±350 c 

735±155 d 

335±150 d 


340±120 c 

380±130 c 

280±160 c 

N e (O III) 

2520±540 e 

4200±700 c 

730±50 d 

140±50 d 

140±70 d 

870±50 c 

730±60 c 

725±50 c 

N e (N II) 


330i^ c 

uotll 6 

37l| dd 


5120±~ 00 d 

oo d 

oo d 

Ne (N II) 






180±80 f 

650i™ of 

420^” f 

Ionic Abundance 

Ratios (xlO -6 ) 








S++/H+ 


4.3±0.7 

12.2±1.6 

14.2±2.1 

20.8±2.7 

9.0±1.4 

14.1±2.2 

12.8±2.3 

0 ++ /H + 

38±5 

71±9 

80±8 

50±6 

82±12 

130±14 

275i31 

227±25 

N ++/ H + 

]4±2 

36±4 

48±5 

27±4 


68±8 

127+14 

117±13 

Ionic Abundance 

Ratios 








N ++/0++ 

0.36±0.04 

0.51+0.02 

0.60±0.04 

0.53i0.07 


0.52±0.03 

0.46±0.02 

0.52±0.02 

o ++ /s ++ 


16.7±2.2 

6.53±0.55 

3.56±0.57 

3.95±0.52 

14.5±1 .2 

19.5±1.9 

17.8±1.8 

N++ /S ++ 


8.47±0.70 

3.90±0.29 

1.89±0.29 


7.54±0.38 

9.04±0.36 

9.18±0.37 

N+/S ++ 


0.63±0.20 

1.90±0.52 

1.75+0.45 

14.3±1.3 

1 .83±0.37 g 

3.98+1. 67 g 

2.30±1.38 6 

N++/N+ 


13.3±4.0 

2.06+0.45 

1.08±0.27 


4.11±0.72 6 

2.27±0.84 s 

3.99±l,46 g 

Column Density Ratios 








0++/S++ 


9.55±0.43 

6.56±0.22 

8.60±0.87 

3.95+0.52 

5.64i0.22 

10.2±0.5 

6.89±0.47 

N++/S++ 


4.85±0.23 

3.92±0.24 

4.56±0.67 


2.94±0.15 

4.72±0.18 

3.56±0.19 

N + /S ++ 


4.62±0.58 

9.96±1 .29 

15.7±2.4 

14.3+1.3 

4.61±0.36 h 

3.29±0.37 h 

1.95±0.50 h 

N++/N+ 


1 05±0.13 

0.39±0.08 

0.29±0.09 


0.64±0.05 h 

1.44±0.14 h 

1.27±0.14 h 

N+/S++ 






3.39±0.41 g 

2.31±0.37 s 

1.55±0.34 5 

N++/N+ 






0.87±0.12 E 

2.04rt0.34 g 

2.30±0.52* 

Adopted Abundance Ratios 








0++/S++ 



6.53±0.55 

3.56±0.57 

3.95±0.52 

14.5±1.2 

19.5±1.9 

17.8±1.8 

N ++/ S ++ 



3.90±0.29 

1.89±0.29 


7.54±0.38 

9.04±0.36 

9.18±0.37 

N+/S++ 


2.63±1.99 i 

5.93±4.03 i 

$.7l±QM 1 

14.3±6. 0 s 

2.96±0.38J 

3.02±0.80 j 

2.59±0.74J 

N ++ /N + 


1.05db0.13 k 

1. 23±0.83 i 

0.69±0.39 i 


1.43± 0.30* 

1.44±0.33 5 

1.89*0.53* 

(N++N ++ )/S ++ 

9.3±3.8 

9.8±4.3 

10.6±7.3 

>14+6 

10.5±0.8 

12.1±1.2 

11.8*1.1 


a Center position, data from Colgan et al. (1993). 

b Radio flux in the aperture used for the FIR line measurement, scaled to the value appropriate for v — 5 GHz. 

c Elect.ron density computed using the 45" aperture data. 

d Electron density computed using the 60" aperture data. 

e Electron density computed using the 80" aperture data. 

f Electron density computed using the 45" aperture data for the 122 /xm line and reducing the 205 /im line flux by the ratio 
of the aperture diameters. 

g Computed using the 45" aperture data for the 122 /xm line and reducing the 205 /xm line flux by the ratio of the aperture 
diameters. 

h Computed using the 60" aperture data for the 122 /xm line and reducing both the 122 /xm and the 205 /xm line fluxes by 
the ratio of the aperture diameters for comparison to the [N III] and [S III] lines. 

Average of column density ratio and ionic abundance ratio; the error is increased for G333.6-N3 where N e { N II) was not 
measured. 

j Average of ionic abundance ratios for all 122 /xm line measurements plus column density for the 45" aperture 122 /xm line 
measurements. 
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Table 3. Estimated Abundances 


Element 

G333.6-0.2 

W43 

Orion Nebula 

Estimated S/H* 

13.4±0.8 x 10“ 6 

12.1±1.1 x 10~ 6 

7.0 x lO" 6 b 

Estimated N/S* 

9.3±2.1 

10.3±0.4 

9.71 b 

Estimated N/H c 

12.5x 10~ 5 

12.5X10' 5 

6.8xl0 -5 d 

Estimated Ne/H 

3.4xl0~ 4 e 

1.1 xl0~ 4 f 

l.OxlO -4 b 

Estimated O/H 

3 to 4xl0~ 4 « 

3 to 8xl<r 4 * 

4.0x 10 -4 d 


a Assumed fractional ionization <S ++ /S>= 0.9. 

b Simpson et al. (1998) revised foT new collisional excitation cross sec- 
tions (see text). 

C N/H = N/S x S/H. 
d Rubm et al. (1991). 

e Estimated from [Ne II] 12.8 //m measurements of Fujiyoshi et al. 
(1998) and radio continuum measurements of Fujiyoshi (1999). 

f SCR£H. 

g O/H = 0*^ + /H + x ICF Mo6el . 
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Table 4. G333.6-0.2 - Comparison of Models to Observations 


Observed Fluxes and Ratios 5 


Model for G333.6-0.2 Core 5,6 Model for North Positions 5 ’ 6 * 


Line Flux or Ratio Center N1 N2 N3 Center N1 N2 N3 N1 N2 N3 


[O III] F 52 fim 

161 c 

92 

[O III] Fgg /rm 

27 c 

39 

[O III] f 52 /f 88 

6.06 

2.37 

[S III] F 19 (it n 

412 c 

270 

[S III] F 33 fim 

158 c 

228 

[S III] F19/F33 

2.61 

1.19 

[N II] F, 22 fim 

2.5 C 

7.4 

[N II] F 2 o5 (im 

0.44 c 

1.97 

[N II] F] 22 /F 2 o5 

5.7 

3.7 

[N III] F 5T (im 

52 c 

55 

[N III] F 57 /[N II] F 122 

21.0 

7.4 

[N III] F s7 /[0 III] F 52 

0.32 

0.60 

(Ne III] F 36 /im c 

17 


[Ne II] F 12.8 Mm d 

1300 


[S IV] F 10.5 (im e 

100 


[Ar III] F 9.0 (im f 

350 


[Ar III] F 9.0 /im 8 

83 


[Ar II] F 7.0 /xm s 

74 


He+/H+ RRLs h 

0 045 


[O III] F 52 (im' 

194 


[O III] Fg 8 fim' 

41 


[N III] F 57 A im 1 

50 



14 

7 

350 

0.0005 

15 

8 

54 

0.0001 

0.96 

0.96 

6.50 

4.12 

61 


869 

273 

79 

47 

351 

243 

0.78 


2.48 

1.13 

2.6 


9.2 

13.7 

1.46 

1.66 

1.5 

2.4 

1.8 


6.2 

5.8 

13 


115 

0.006 

5.0 


12.6 

0.0004 

0.89 


0.33 

12 



8.5 




1793 




95 




251 




150 




13 




0.044 




359 




56 




144 



0 

0 

117 

1 

0 

0 

0 

67 

1 

0 



1.75 

1.61 


50 

14 

154 

52 

16 

72 

27 

192 

85 

33 

0.70 

0.52 

0.80 

0.61 

0.48 

5.6 

2.8 

6.8 

6.9 

3.5 

1.4 

0.9 

1.9 

2.0 

1.4 

4.1 

3.0 

3.7 

3.4 

2.5 

0 

0 

62 

1.9 

0 

0 

0 

9.1 

0.28 

0 



0.53 

1.49 



& Fluxes xlO 18 W cm 2 , corrected for extinction. 

b Model Abundances (x!0“ 6 ): C/H = 400, N/H = 125, O/H = 500, Ne/H = 340, S/H = 13.4, Ar/H = 4.5. 

c 0.75 / Aperture, Colgan et al. (1993). 

d 0.75 ; Aperture, Geballe et al. (1981). 

e 0.50' Aperture, Geballe et al. (1981). 

f 0.67' Aperture, Geballe et al. (1981). 

s 0.35' Aperture, Cohen et al. (1989), Simpson et al. (1995b). 

h Radio recombination lines measured in 2.3' aperture, McGee Newton (1981). 

‘80" Aperture for ISO LWS. 



Table 5. W43 - Comparison of Models to Observations 


Ratio Observed Ratios Model D a Model S a Model T a Comments 


He + /H + RRLs 

0.076 b 

0.075 

0.073 

[S III] 19/33 m™ 

0 78, 0.82, 0.72 

1.04 

1.03 

[O III] 52/88 M® 

2.6, 2.3, 2 3 

2.30 

2.33 

[N II] 122/205 Mm 

4.3, 7.1, 6.2 

6.29 

6.28 

S ++ /S 


0.919 

0.953 

Q++/S++ « 

16.33±0.89 

0.583 

0.391 

N ++ /0 ++ 

O.oOiO.Ol 

1.18 

1.58 

N ++ /N + 

1.5 to 3.0 

1.72 

1.45 

N+/S++ 

2.90±0.31 

0.400 

0.428 

N++/S++ 

8.61=^0.21 

0.686 

0.620 


0.075 Good r eff 

1 .02 Model N € too high 

2.33 Good N e 

6.19 Good N e 

0.965 Implies S/H = 11.8, 11.4, or 11.3xl0~ 6 
0.221 Implies O/H = 3.3, 4.8, or 8.3xl0 -4 d 
3.05 Implies N/O = 0.42, 0.32, or 0.16 e 
1.86 Compare uniform He + /H + 

0.363 Implies N/S = 7.25, 6.78, or 7.99 f 
0.673 Implies N/S = 12.55, 13.89, or 12.79 f 


a The stellar atmosphere used for the Model D is the sum of 0.8 times the D-35 model atmosphere from Pauldrach 
et al. (2001) plus 0.2 times the D-40 model atmosphere, the stellar atmosphere used for Model S is the sum of 
0.6 times the S-35 model atmosphere plus 0.4 times the S-40 model atmosphere, where D-35, D-40, S-35, and S-40 
stand for Dwarf and Supergiant models with T& equal to 35,000 K and 40,000 K, respectively, and the stellar 
atmosphere used for Model T is the TLUSTY model atmosphere from Lanz k Hubeny (2003) for T e $ = 35,000 K 
and log g = 3.50. The model abundance ratios (xl0~ 6 ) are C/H — 400, N/H = 125, O/H = 500, Ne/H = 110, S/H 
= 12.1, and Ar/H = 4.5. 

b Weighted average of measurements of Churchwell et al. (1978), Lichten, Rodriguez, & Chaisson (1979), Lockman 
& Brown (1982), Peimbert et al. (1992). 

c Model 0 ++ /S' f+ is the density weighted ratio of <0 +_( "/0>/<S ++ /S>, etc. (SCREH). It should be multiplied 
by the relevant abundance ratio for comparison with the observed ratios. 

d O/H = O/S x S/H. 

e N/H = 1.4, 1.5, or 1.4 xlO" 4 if O/H = 3.3, 4.8, or 8.3 xl0“ 4 . 

f The average N/S = 9.9, 10.3, or 10.4, N/H = 1.17, 1.18, or 1.17 x!0~ 4 , and O/H = 2.8, 3.7, or 7.2 xl0“ 4 . 


